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Abstract 

A new thermodynamic phase resulting from the competition between a com- 
mensurate potential and disorder in interacting fermionic or bosonic systems 
is predicted. It requires interactions of finite extent. This phase, intermediate 
between the Mott insulator and the Anderson insulator, is both incompress- 
ible and has no gap in the conductivity. The corresponding phase is also 
predicted for commensurate classical elastic systems in presence of correlated 
disorder. 
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The interplay between disorder and interactions gives rise to fascinating problems in 
condensed matter physics. Although the physics of disordered noninteracting systems is 
by now well understood, when interactions are included the problem is largely open, with 
solutions existing for one dimensional systems or through approximate scaling |3|-[5| or 
mean-field methods |6j . A particularly interesting situation occurs when the non-disordered 
system possesses a gap. This arises in a large number of systems such as disordered Mott 
insulators ]7|-|TU|], systems with external (Peierls or spin-Peierls systems [II|) or internal 



commensurate potential (ladders or spin ladders [^2 -[15], disordered spin 1 chains 
This situation is also relevant for classical problems such as elastic systems subjected to 
both periodic potential and correlated disorder, as encountered e.g. in vortex lattices in 
superconductors fl8fl . 

Although in some cases an infinitesimal disorder can suppress the gap |7j due to Imry-Ma 
effects, in most finite amount of disorder is needed to induce gap closure. In the latter 

case, the complete description of the gap closure and of the physics of the resulting phases is 
extremely difficult with the usual analytic techniques such as perturbative renormalization 
group (RG), due to the absence of a weak coupling fixed point. In d — 1, where attempts 
in solving this problem could be made, it was believed |8||| that two phases existed: a weak 
disorder phase where the gap is robust and the system has all the characteristics of the pure 
gaped insulator, and a strong disorder phase where the gap is totally washed out by disorder 
and the system is a simple compressible Anderson insulator. However the techniques used 
so far are either restricted to special points || or suffer from serious limitation: the simple 
perturbative RG of [@] has to be used outside its regime of validity to describe the transition 
between two strong coupling phases. Furthermore up to now mostly onsite interactions have 
been studied. 

In this Letter we thus reexamine this problem using better suited methods that capture 
some nonperturbative effects: a variational calculation and a functional renormalization 
group method. We focus first, for simplicity, on one dimensional interacting spinless fermions 
in the presence of a commensurate potential |19[ . As argued below we expect similar physics 



to hold for systems with spins. Our main finding is that, in addition to the above mentioned 
phases, an intermediate phase exists, as shown in Fig. [I]. 

Quite remarkably although this phase possesses some of the characteristics of an Ander- 
son insulator, namely a non-zero ac conductivity at low frequency, it remains incompressible. 
Since it is also dominated by disorder, we call it a Mott glass. We discuss here the physical 
characteristics of this novel phase and its interpretation for fermionic as well as bosonic 
systems. We argue that such a phase is not specific to d = 1 and discuss its generalization 
to higher dimension both for quantum and classical systems. 

Disordered interacting spinless fermions submitted to a commensurate periodic potential 
are described in one dimension by 



H = dx 



where ± denote fermions with momentum close to ±A;f, as is standard in d = 1. V is the 
interaction, /j,(x) an on site random potential and g the strength of the periodic potential 
opening the gap. In d = 1 it is convenient to use a boson representation of the fermion 
operators pi. This leads to the action 
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where a is a short distance cutoff (of the order of a lattice spacing) and the field is related 
to the total density of fermions by p(x) = — V4>/tt. The interactions are totally absorbed in 
the Luttinger liquid coefficients v (the renormalized Fermi velocity) and K, a number that 
controls the decay of the correlation functions. The noninteracting case (V = 0) gives K = 1 
and v = Dp. We focus in the following on the repulsive (V > 0) case for which < K < 1 
depending on the strength and range of V. For weak disorder one can separate in the random 
potential the Fourier components close to q ~ (forward scattering) and q ~ 2kp (backward 
scattering). We have retained here, for simplicity, only the backward scattering. Forward 



scattering, studied in [21], does not lead to qualitative changes in the physics presented here. 



The backward scattering is modeled by the complex gaussian random variable £ such that 
= W5(x-x'). 

Perturbatively both the disorder (in the absence of commensurate potential) and the 
commensurate potential (in the absence of disorder) are relevant (respectively for K < 3/2 
and K < 2). Let us now study this problem using a variational method |22j. We first 
use replicas to average over disorder in ([!]). We then approximate the replicated action by 
the best quadratic action So = \ L u ^q^G^b (q, u} n )(j)g tU}n where the G are the variational 
parameters and a, b replica indices . Skipping the technical details which can be found 
in I^TJ, we find that the saddle point solution yields 

vG-\q,u n ) = ^=(u 2 n + v 2 q 2 ) + m 2 + E^l - 6 nfi ) + I{oJ n ) (2) 

where G c = J2b G a b is the connected Green function. The parameters m, Si and the func- 
tion I{oj n ) obey closed self consistent equations. The important physical quantities are 
simply given in term of (|2|): the static compressibility reads k = lim^o lim^o q 2 G c (q, uj) 
while the conductivity is given by the analytical continuation to real frequencies: cr(u) = 
[uj n linig^o G c (q, c^n)]!^ Note that both quantities are stemming from the same prop- 

agator but with different limits q — > and uj — > 0. 

Since we expect the physics to be continuous for small enough K (i.e. repulsive enough 
interactions), one can gain considerable insight by considering the classical limit h — > 0, 
K — > keeping K = K/h fixed. In this limit one can solve analytically the saddle point 
equations and compute m, Si and I(ui n ). The resulting phase diagram is parameterized 
with two physical lengths (for K — > 0): The correlation length (or soliton size) of the 
pure gapped phase d = ((AgK)/(av )) -1 ^ 2 0] and the localization (or pinning) length l = 
((av ) 2 /(l6WK )) x / 3 |I|J25| in the absence of commensurability. We find three phases as 
shown in Figure [I]: 

Mott insulator [MI]: At weak disorder we find a replica symmetric solution with Si = 
but with m ^ (m depends on the disorder). m^0 leads to zero compressibility k = 0. m 
defines the correlation length £, with £ 2 = v 2 /(nKm 2 ), in the presence of both the disorder 
and the commensurate potential. The effect of disorder is to increase £ compared to the pure 
case, since it reduces the gap created by the commensurate potential. We find that £ is given 



by (d/O 2 exp[±(£// ) 3 ] = 1. We obtain a(cu) = if uj < u c = m^Jl + A - 3(A/2) 2 /3 where 
A = (£//q) 3 - The physics of this phase is similar to the simple Mott insulator. However the 
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gap in the conductivity decreases, when disorder increases, and closes for A = 2. For A > 2 
(d/lo > 0.98) the RS solution becomes unphysical even though the mass m remains finite at 
this transition point. 

For stronger disorder one must break replica symmetry. In the absence of commensurate 



potential such a solution describes well the d — 1 Anderson insulator p3|| , in which Si ^ 0. 
Here, however two possibilities arise depending on whether the saddle point allows for m/0 
or not: 

Anderson Glass [AG]: For large disorder compared to the commensurate potential d/l > 
1.58, m = is the only saddle point solution. In this case one recovers exactly the solution 
of the Anderson insulator with interactions but no commensurate potential. We call it 
Anderson glass to emphasize that this phase is dominated by disorder (the corresponding 
phase for bosons, also described by ([I]), is the Bose glass [1]J5[])- Within the variational 
approach, the AG has a finite compressibility (identical to the one of the pure system 
k = ttK/v) and the conductivity starts as o~{uj) ~ uj 2 showing no gap. Physically this is 
what is naively expected if the disorder washes out completely the commensurate potential. 
While the MI and AG were the only two phases accessible by previous techniques we 
find that an intermediate phase exists between them. 

Mott Glass [MG]: For intermediate disorder 0.98 < d/l < 1.58 a phase with both S x ^ 
and m exists, m and £i define two characteristic lengths in the intermediate phase, 
m 2 + £i remaining constant in the MI and MG [^TJ. On the other hand m varies and vanishes 
(discontinuously within the variational method) at the transition from MG to AG. The MG 
is thus neither a Mott nor an Anderson insulator. In particular, the optical conductivity has 
no gap for small frequencies o~{uj) ~ u 2 , while due to m ^ 0, the system is incompressible 
k = 0. These properties are shown in Fig. ||. 

This result is quite remarkable since by analogy with noninteracting electrons one is 
tempted to associate a zero compressibility to the absence of available states at the Fermi 
level and hence to a gap in the conductivity as well. Our solution shows this is not the 
case, when interactions are turned on the excitations that consists in adding one particle 
(the important ones for the compressibility) become quite different from the particle hole 
excitations that dominate the conductivity. 

Physical arguments are also in favor of the existence of the Mott Glass, both for systems 
with or without spins. Let us consider the atomic limit, where the hopping is zero. If 
the repulsion extends over at least one interparticle distance, leading to small values of K, 
particle hole excitations are lowered in energy by excitonic effects. For example for fermions 
with spins with both an onsite U and a nearest neighbor V the gap to add one particle is 
A = U/2. On the other hand the minimal particle-hole excitations would be to have the 
particle and hole on neighboring sites (excitons) and cost A p .h. = U — V. When disorder is 
added the gaps decrease respectively |26| as A — > A — W and A p .h. — > A p .h. — 2W. Thus the 
conductivity gap closes, the compressibility remaining zero |[27|| . According to this physical 
picture of the MG, the low frequency behavior of conductivity is dominated by excitons 
(involving neighboring sites). This is at variance from the AG where the particle and the 
hole are created on distant sites. This may have consequences on the precise low frequency 
form of the conductivity such as logarithmic corrections. When hopping is restored, we 
expect the excitons to dissociate and the MG to disappear above a critical value K > K*. 
Since finite range is needed for the interactions, in all cases (fermions or bosons) K* < 1. 
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In addition we expect K* < 1/2 for fermions with spins. Similar excitonic arguments 
should also hold in in two- or three-dimensional bosonic and fermionic Mott insulators 
provided some finite range of interaction is taken into account. In higher dimension, since 
disorder has a weaker impact on the transport properties, one expects that the important 
change in the conductivity occurs at the transition between the MI and MG, whereas the 
compressibility would become non zero only for stronger disorder (transition MG to AG). 
Numerical investigations would prove valuable. Small gaps facilitate the observation of MG 
physics (see Fig. |]), making the study of systems already close to a metal - Mott insulator 
transition particularly interesting . 

The properties that we have obtained are thus quite general, depending only on the two 
gaps closing separately. Given the mapping between a d dimensional quantum problem and 
a d + 1 classical one our study also applies to commensurate classical systems in presence of 
disorder correlated in at least one direction (here the imaginary time r). (0) can be general- 
ized to any dimension to describe a classical elastic system where becomes a displacement 
field u. It applies to systems with internal periodicity, such as a crystal or a charge density 
wave (with 2(f) = Kqu for reciprocal lattice vector Kq) or to non periodic systems such as 
interfaces. The periodic potential making the system flat while the disorder makes it rough. 
For these systems a functional RG procedure (FRG) in d = 4 — e (i.e near 4 + 1-dim systems) 
can be used. For uncorrelated disorder in the absence of an external periodic potential an 
internally periodic system is described ino?<4byaT = "Bragg glass" fixed point |p8| 



Adding a periodic potential cos(p0) as a perturbation, a transition was found at T = [29 
between the Bragg glass (for p > p c (d)) (periodic potential irrelevant) and a commensurate 
phase (for p < p c (d)) (periodic potential relevant). Such a T = transition also exists for 
correlated disorder 

To confirm the existence of the MG phase it is necessary to study the phase where 
the periodic potential is relevant (since p = 2 < p c (d) here). Since this goes beyond the 
perturbative FRG analysis, we consider the toy model where the cosine is replaced by a 
quadratic term (a good approximation when the periodic potential is relevant) defined by 
the energy: 

Y = tJ ddxdT ll(<Vu) 2 + c u (d T u) 2 + m 2 u 2 ) + V(x,u(x,t))] (3) 

where r is the coordinate along which disorder is correlated (e.g. the magnetic field for 
vortices with columnar defects), T is the classical temperature (h for the quantum problem) 
and the gaussian disorder has a correlator V(x, u)V(x', u') = 5 d (x — x')R{u — v!) (for a 
periodic problem, R(u) is itself periodic), c and c 44 are the elastic moduli, analogous to 1/K 
for the quantum problem ([!]). This d + 1 dimensional model can be studied perturbatively 
(for small m and disorder). For m = and T = 0, both for uncorrelated |3(J and correlated 



disorder a cusp-like nonanalyticity in the renormalized disorder R(u) develops beyond 
the Larkin length R c (corresponding to metastability and barriers in the dynamics) |j32 



Remarkably, we find that this feature persists even when m > 0, while one usually expects 
that a mass smoothes out singularities. This can be seen from the RG equation for A(u) = 
-R"(u): 

^A(n) = eA(u) + T}A» + /,(A"(«)(A(0) - A(n)) - A'(u) 2 ) (4) 
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with fi = ^2(1 + fie 2l )~ 2 , A* = m2fl2 and at zero temperature TJ = 0. Integrating the closed 
equation for A"(0) one finds that the cusp persists (i.e — A" =+oo (0) = +oo) provided that 
R c < R*(fM) (with ~ l/y/U for small //) while it is washed out (-Af =+oo (0) < +oo) 

for weaker disorder. Our FRG study shows that this T = transition in the renormalized 
disorder, exists both for correlated and uncorrelated disorder. For uncorrelated disorder 
our findings are of interest for the question of the existence of an intermediate "glassy flat 
phase". In that case, however, no sharp signature of this transition exists in two point 
correlation functions and a physical order parameter remains to be found, which makes 
the existence of such a phase still controversial On the contrary for correlated 

disorder the transition seen in the FRG has much stronger physical consequences. Because 
of the lack of rotational invariance (in (x,t)) the existence of the cusp and the transition 
directly affects two point correlation functions. Indeed, the tilt modulus c 44 renormalizes 
as <9;C44 = — /(/)A"(0)c44. Integrating at T = one finds that C44(/ = +oo) is finite for 
R c > R* c {n) but that it is infinite for R c < R* c (n). Furthermore we find (see [|2l] for details) 
that for correlated disorder a small temperature (T\ > 0) does not affect the transition (the 
phase where c 44 (Z = +oo) = oo survives) whereas for uncorrelated disorder the effective 
temperature goes to a constant (~ /i 1 " 6 / 2 ) washing out the cusp. Thus the toy model 
for correlated disorder exhibits at low temperature, within the FRG, a transition between 
two phases. The first one is identified with the MI, where the mass (commensurability) 
destroys the metastability and restores isotropy in x, r at large scale. The second one is 
the MG, which is glassy with metastable states despite the presence of the mass. It shares 
some properties with the AG such as C44 = 00. This implies non analyticity of the Green's 
functions in frequency. The AG itself corresponds to the phase where the periodic potential 
is irrelevant (m = 0). This provides strong evidence for the intermediate phase proposed in 
this paper, which, besides electrons, should be obtained in classical commensurate elastic 
systems with disorder. 
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FIGURES 




FIG. 1. The three phases as defined in the text in the disorder W, commensurate potential 
g plane for a fixed value of the interactions. Both MI and MG are incompressible. Both MG and 
AG have no gap in the conductivity. 
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FIG. 2. Conductivity in the MI (solid line) and MG phases (dashed line). Insert: density n 

versus the chemical potential fi. 
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